Is transport in time-dependent random potentials universal ? 
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The growth of the average kinetic energy of classical particles is studied for potentials that are 
random both in space and time. Such potentials are relevant for recent experiments in optics and 
in atom optics. It is found that for small velocities uniform acceleration takes place, and at a later 
stage fluctuations of the potential are encountered, resulting in a regime of anomalous diffusion. 
This regime was studied in the framework of the Fokker-Planck approximation. The diffusion 
coefficient in velocity was expressed in terms of the average power spectral density, which is the 
Fourier transform of the potential correlation function. This enabled to establish a scaling form 
for the Fokker-Planck equation and to compute the large and small velocity limits of the diffusion 
coefficient. A classification of the random potentials into universality classes, characterized by the 
form of the diffusion coefficient in the limit of large and small velocity, was performed. It was 
shown that one dimensional systems exhibit a large variety of novel universality classes, contrary to 
systems in higher dimensions, where only one universality class is possible. The relation to Chirikov 
resonances, that are central in the theory of Chaos, was demonstrated. The general theory was 
applied and numerically tested for specific physically relevant examples. 



Dynamics in potentials which are random both in space 
and time were subject of many sophisticated studies for 
nearly 100 years [iRjl- The response to forces resulting 
of such potentials typically differs from ordinary diffu- 
sion. Specifically, the diffusion coefficients predicted by 
such mechanisms sensitively depend on the velocity of 
the particles. Also, if the potential is time-dependent, 
the energies of the particles will not be constant. The 
existence of such 'anomalous diffusion' has been demon- 
strated for classical dynamics with spatially and tempo- 
rally fluctuating potentials [5|-[lOJ. These works claim 
universal behavior in the sense that for generic random 
potentials the diffusion coefficient exhibits a universal 
power-law dependence on instantaneous velocity w, such 
that D{v) ^ |w|~'^ as \v\ -^ oo. This in turn implies 
that asymptotically in time the average velocity satisfies 
(w^) '^ t^/^ . In the present Letter we will demonstrate 
that this picture should be revised, and we will intro- 
duce such a revision. The potentials we study here are 
continuous in space and are fundamentally different from 
lattice models where the velocity (momentum) is inher- 
ently bounded. 

In the present work we consider the classical dynam- 
ics of a particle in potentials that are random both in 
space and in time, emphasizing the spreading of the ve- 
locity acquired by the particle, as time evolves. This is 
a fundamental problem, which was motivated by experi- 
ments in optics [11] and in atom optics 12ill3|, , where 
the random potential is introduced by transforming an 
intensity pattern into an effecti ve p otential, for the light 
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Such potentials are 
naturally described in terms of the Fourier spectrum of 
the waves inducing them, and their spectral coefficients 
are assumed to be independent random variables. In ad- 
dition to the fundamental interest, the calculations of 
the present work are also relevant for spreading of waves 
in random time-dependent potentials, because at least 



in the regime of large velocities (short wave-length) it is 
generally believed that a classical picture is appropriate. 
Therefore, we may expect that the classical results, which 
are presented in this Letter are relevant also for the wave 
experiments (e.g., [Ill4l3j ) . We will consider a general 
stationary (both in space and time) random potential, 
which is also isotropic on average. Such potentials are 
conveniently described by their Fourier components, 

y(i)(x,t)= / y(k,w)expi(k-x-wt)dkda;-hc.c, (1) 



where V (k, w) is a random field chosen such that 
the distribution of T^^^^ (x, t) is stationary both in 

time and space, in such case (V {\i,(jj)V {\i! ,(jj')) = 

VJj^iJ (k — k') (5 (cj — cj'), where (.) denotes the ensemble 
average. For time dependent potentials the velocity of 
the particles may grow, since energy is not conserved. 
The dominant mechanism of the growth of the velocity 
is via Chirikov resonances |16j . namely, resonances be- 
tween the particle dynamics and the external driving. 
Those resonances occur when the phases in ([T]) are sta- 
tionary. For a given velocity v = x, this happens for 



k • V = 



(2) 



In one dimensional systems this condition reduces to 
kv — id, namely, when the particle just 'surfs' on one 
of the waves composing ([1]) . The growth of the velocity 
occurs only when the particle 'jumps' to a wave travel- 
ing with a nearby and larger velocity than the velocity of 
its initial carrier. This provides an intuitive picture for 
the understanding of transport in phase-space, since the 
growth of the velocity is limited to regions in phase-space 
where the density of Chirikov resonances is non-zero. In 
dimensions two and higher particles may increase their 
velocity beyond the velocity of the wave by traveling in 
a direction that is not parallel or even perpendicular to 



the wave propagation direction. This essential difference 
wiU yield distinct much richer transport behavior in one 
dimensional systems than in higher dimensions, as will 
be explained in what follows |17| . 

Since, in this work we consider only potentials contin- 
uous in both time and space it is natural to define the 
characteristic length and time scales of the potential, l^ 
and It, correspondingly. Those, scales are defined such 
that the variation of the potential over them is limited, 
and the particle will encounter an almost constant force. 
Note, that It and Ix are generally different from the ex- 
ponential decay rates of the correlation function of the 
potential. If the force is weak in the sense, F < Ix/lt, 
particles with velocities ■;; < Ix/h (as will be assumed in 
the present work) will experience an uniform acceleration 
up to time t < It- For longer time-scales the variations of 
the potential become apparent for the particle, and the 
velocity will exhibit anomalous diffusion. 

To calculate the anomalous diffusion coefficient we 
invoke the Fokker-Planck approximation, assuming 
that the force is sufficiently weak and the decay 
of the potential correlations, C (xi — X2, ti — t2) = 
{V {xi,ti)V {x2,t2)) , is sufficiently rapid, so that the 
velocity can be considered constant on the time scale 
where the correlation function is appreciable. For sta- 
tionary random potentials it is convenient to use the av- 
erage power spectral density (PSD), S{k,uj), defined as 
the Fourier transform of the correlation function (Weiner- 
Khinchin theorem). 



C (x, t)= duj dkS (k, uj) exp i (k • x - cut) 



(3) 



The benefit of this representation is twofold. First, in 
experiments Ill4l3lll8l.ll9j the PSD is the naturally con- 
trolled rather than the correlation function. Second, the 
analytical relations obtained in this representation are 
more transparent. For a potential which is not only sta- 
tionary, but also has an isotropic PSD, and initial dis- 
tribution which is isotropic as well, the Fokker-Planck 
equation for the velocity is effectively one dimensional 



f=(„-<-«iL„-.z,(„)|)p. (4) 

where P (v, t) is the probability density and D (v) is the 
diffusion coefficient given by 



D{v) 



dk (k-vyS(k,k-w) 



(5) 



Note, that if the PSD has some typical scales, /cq and wq 
then the Fokker-Planck equation is invariant under the 
transformation of variables, 



t' ( nV„'^ 

LOO 



which allows the rescaling of the Fokker-Planck equa- 
tion to a universal form (where Vq is the amplitude of 
the potential). For dimensions two and higher, changing 
variables to y = v cos 9 and expanding the angular part 
of the the integrand in ((S]) gives the asymptotic behavior. 



Div)^^, 



with 



D3 = 2Sd dy dky^k''+^S{k,ky), 

Jo Jo 



(7) 



(8) 



where Sd is the surface of a d— dimensional hyper-sphere. 
Using this asymptotics an asymptotic scaling solution to 
the Fokker-Planck equation may be obtained, P{v,t) = 
i~''/'^g (w^/t), which yields the growth of the mean ki- 
netic energy as, i (w^) '~ t'^'^. This behavior is consid- 
ered in the literature as universal [5|, [Zl, [sj, |lO| for any di- 
mension given that the correlation function of the poten- 
tial is sufficiently differentiable. For dimensions two and 
higher as shown above this is indeed the case, however 
for one dimensional systems other behaviors are possible. 
For one dimensional systems new possibilities arise 
from the fact that the angular part in ([5|) is missing and 
therefore the diffusion coefficient. 



Div) 



k^S{k,kv)dk. 



(9) 



(6) 



may acquire various asymptotic behaviors dictated by 
S (fc, w). Consider a diffusion coefficient which decreases 
with velocity faster than any power law, then clearly the 
asymptotic expansion in w~^ is not useful, since it will 
produce a nil result. A simple example is given by a 
PSD with the property, S{k,uj) = uj/k > timax, 
which renders the diffusion coefficient zero for large ve- 
locities, D (v) = V > Wmax, as demonstrated in [17| 
for a specific example. We will now demonstrate this be- 
havior, using an experimentally relevant potential which 
is proportional to the intensity, V {x,t) — \U {x,t)\ 
of some complex field, which in turn is a superposi- 
tion of waves U (x, t) — f dkU (k) exp i {kx — uj (fc) t) , 
with some dispersion relation u;{k). These potentials 
appear in experiments with neutral atoms [1^ ll3| and 
in some experiments in optics [ll|, [l8|, [l9[, where the 
variation of the refractive index (which plays the role 
of the potential) of a photosensitive material is propor- 
tional to the intensity of light [ij, [l^- For simplicity 
we will assume the dispersion relation to be uj (fc) = 
fc^/2 (which naturally appears in the framework of the 
paraxial approximation in optics), and / (fc) to be the 
probability density of the wave numbers k. For this 
type of potentials one can readily obtain both the PSD, 
S{k,uj) = Vi \k\-^ / (I _ I) / (I + I) ,which is singu- 
lar for k 



0/191/(9 



0, and the diffusion coefficient, D {v) 
v) f {q — v) dq. It is clear that the decay 
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Figure 1: A log-log plot of average squared velocity as a 
function of time for one dimensional system with a poten- 
tial V (x,t) = \U {x,t)\ and a uniform distribution of wave- 
numbers. The blue dots represent the result of Monte-Carlo 
simulation the corresponding Newton equations averaged over 
20 realizations and the black solid line is the numerical solu- 
tion of the Fokker-Planck equation for the velocity Q. The 
dashed black and red lines are guides for the eye with the 
corresponding slopes of 2 , 1 and 0. The initial condition 
was a narrow distribution of velocities for the Fokker-Planck 
Eq. and x = v = Q for the Monte-Carlo calculation. The 
parameters used for this simulation are, Vb = 10"*, ku — 0.1. 



of the diffusion coefficient with the velocity is dictated by 
the decay of the probability density of the wave-numbers, 
/ (fc). An explicit expression may be obtained for exam- 
ple for a uniform distribution of wave-numbers (in the 
interval [-kR^kn]) giving [13, 



D{v) = \^ ^ ' 



\v\y 



\v\ < kR 

\V\ > kR 



(10) 



The resulting dynamics is demonstrated in Fig. [TJ The 
regime of a unit slope corresponds to regular diffusion 
and the asymptotic regime corresponds to an absence of 
diffusion in velocity. 

Another explicit expression may be obtained for a 
Gaussian distribution of wave-numbers, giving a diffu- 
sion coefficient 



D{v) = 2Vo^e^p{-vykl). 



fll) 



The predictions of this equation are compared with 
Monte-Carlo simulations and the results are presented 
in Fig. [5] The correlation function for this potential can 
be also obtained using ^, 



Cix,t)^Vo^ 



1 



VTTW 



exp 



l + fc4i2 



(12) 



Note, that the correlation function ([T2|) is infinitely dif- 
ferentiable, decaying fast in position and slowly decaying 
in time. 

These examples demonstrate asymptotic behavior 
which is very different from ([7]) described in previous 



Figure 2: Same as Fig. [T] but with a Gaussian distribution 
of wave-numbers. The dashed black lines are guides for the 
eye with the corresponding slopes of 2 and 1. The parameters 
used for this simulation are, Vb = 10~^, ka = 0.1. 



studies [5|, |8|, |lO| . Therefore new universality classes of 
potentials should be defined. These universality classes 
are characterized by a growth in velocity which is sub- 
diffusive such that the moments of the velocity distribu- 
tion grow slower than any power law. It is important 
to note that these new universality classes are neither 
limited to potentials with a dispersion relation nor to 
correlations that are slowly decaying. For example, the 
PSD, 

S (fc, w) - exp i-ko/kf exp {~uj^/uJq ~ k^/k^) , (13) 

gives rise to a rapidly decaying correlation func- 
tion and a diffusion coefficient which is D (v) ^ 



is some polyno- 



g (v) exp ( - Y 1 + (v/vaf j , where g (u) i 

mial of V. The classification of potentials into different 
universality classes, which differ from ([7]), may be ob- 
tained for a PSD which is not singular and differentiable 
on the line k = 0. This turns out very useful since in 
some experiments the PSD could be precisely controlled, 
which suggest a way to test the various claims of this Let- 
ter. The diffusion coefficient can be expanded in powers 
oi v~^, 



D{v) 






D. 



+ 



(14) 



where 



Dn = 



(n-3)! 



dojw"^ -^zTz — -{k,uj)\k=() n>3. 



dk 



n— 3 



(15) 

If the first Umax derivatives, d^S/dk", vanish on the line 
fc = 0, where rimax is the first derivative that is dif- 
ferent from zero, then the resulting asymptotic behav- 
ior is -D (w) ~ _D„^_^^+3/w""""'+'^. The case where all 
the derivatives vanish is possible only if S'(fc,a;) is non- 
analytic on the line fc = 0. The non-analytic behavior 



may be found either for a function, which is strictly zero 
on some finite strip around the hue fc = 0, or due to 
an essential singularity of S {k,uj) on this line. The first 
case will lead to a diffusion coefficient that will vanish 
for large velocities, for example PH)) (where the support 
of / (fc) is finite). The second case results in a decay 
faster than any power law (jlip. By controlling the an- 
alytical behavior of the PSD on the line fc = 0, one can 
vary the asymptotic behavior of the velocity dependence 
of D (v). The range of variation is from D (v) r^ v~^ , 
through D (v) ^ v~" (with n > 3) to sub-exponential, 
exponential, super-exponential and up-to D (v) — (for 

In this Letter the diffusion coefBcient for the velocity 
is presented in terms of the average power spectral den- 
sity (PSD) (see Eq. ([5])) for stationary potentials that 
are random both in space and time. This representation 
is very natural for stationary (both in time and space) 
potentials, in particular, for potentials which are a su- 
perposition of waves, as they appear in optics and atom 
optics. The simplicity of the expression enabled to ex- 
plore the properties of the diffusion coefficient, to estab- 
lish a scaling form of the Fokker-Planck equation, and to 
discover new universality classes. In particular, we were 
able to calculate explicitly the diffusion coefficient for 
representative examples, relevant for applications both 
in optics and atom optics (Gaussian and uniform distri- 
butions of wave- vectors). We have shown that for di- 
mensions larger than one only one universality class is 
possible in the framework of the Fokkcr-Planck approxi- 



mation, D (v) 



However for one dimensional sys- 



tems new possibilities for large velocity asymptotics were 
also found. It was demonstrated that diffusion in phase- 
space takes place only where Chirikov resonances ([2]) are 
found. 

The main result of this Letter is the classification of 
universality classes in one dimensional systems. In the 
past it was found that in the large velocity limit the dif- 
fusion coefficient depends on the velocity as, D (v) ~ v^'^ 
p, M, la, [lOJ . In the present work we have shown that this 
is always the case for dimensions larger than one, and 
explained the mechanism of this behavior (see discussion 
before (O). However, for one dimensional systems this 
is only one of the possibilities. Generally, the possible 
asymptotic expansion of D (v) in powers of v^^ is: 

1. The first term of the asymptotic expansion of D (v) 
is non-zero, D (v) ^ v^^ . 

2. The first non- vanishing term in the asymptotic ex- 
pansion of D (v) is n then, D (v) ^ v^". 

3. All terms in the asymptotic expansion oi D (v) are 
zero, D (v) < v~", for any a > 0. In particular, 
the diffusion coefficient may be zero for v > Wmax 
(e.g. PU]) ) or non-zero but decreasing faster than 
any power law (e.g., (ITT|) '). 



All these possibilities can be realized in experiments 
with good control over the PSD or more precisely the 
probability density of the wave- numbers, which comprise 
thepotential, /(fc). Unlike statements of other studies 
[a la, IM |lO| , the differentiability of the correlation func- 
tion of the potential is not related to the classification 
into universality classes. Additionally, it was demon- 
strated that the new universality classes do not depend 
on the range of the correlation function. 

For small velocities the diffusion coefficient, obtained 
in the framework of the Fokker-Planck approximation, 
will be generally different from zero and therefore a 
regime of regular diffusion is expected. However, we 
have shown that initially particles will experience a uni- 
form acceleration, and therefore at least for short times 
the Fokker-Planck approximation is invalid, as could be 
clearly seen from all the figures. We have also found 
potentials, which produce diffusion coefficients grow ing 
with velocity (from zero and up to some value) [l7[ for 
which the validity of the Fokker-Planck approximation is 
not satisfied, since the condition that the velocity is con- 
stant during the correlation time is violated. On a longer 
time-scale, the velocity will follow a diffusion equation 
and eventually, it will reach the asymptotic long-time be- 
havior of anomalous diffusion as predicted by the Fokker- 
Planck approximation. 

In the present work the spreading of the velocity distri- 
bution is studied in the framework of the Fokker-Planck 
equation, therefore an obvious question to study is what 
happens when the Fokker-Planck approximation fails. 
Another important issue, which was not addressed in this 
work, is how to obtain analytically the spreading in posi- 
tion. Since in the experiments that have motivated this 
work 



11 



19|, the relevant dynamics is of waves, rather 
than particles, an obvious question to explore is the cor- 
respondence between the classical and wave dynamics. 
This correspondence is shown in Fig. |3l where in addi- 
tion to Monte-Carlo and Fokker-Planck calculations of 
the average squared velocity for particles a numerical 
simulation of a corresponding wave system is presented. 
As is expected there is a reasonable correspondence for 
large velocities. It will be explored in further studies and 
the possibility for its violation is of great interest. 

This work was motivated by the experimental work 
of Liad Levi and Mordechai Segev, whom we thank for 
many stimulating discussions and for providing crucial 
insight for this problem. It is our great pleasure to thank 
Tom Spencer for introducing us to a, [a and to Michael 
Wilkinson for introducing us to [8|, |lO[. Many of the 
results of the present work originated from fruitful dis- 
cussions with Michael Wilkinson during his visit to the 
Technion. The work was supported in part by the US- 
Israel Binational Science Foundation (BSF). 
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Figure 3: Same as Fig. [5] but for a two dimensional system, 
and including the wave dynamics . The dashed black lines 
are guides for the eye with the corresponding slopes of 2 and 
2/5. The dashed (light) red line is the result of a numerical 
simulation of the Schrdinger equation with the same potential, 
averaged over 100 realizations. The initial distributions of 
particle velocities are the same as the initial wavefunction, a 
narrow Gaussian around the origin. The parameters used for 
this simulation are, Vb = 2.06 x 10~*, kn = 0.029 (same as 
experimental parameters in [131 )■ 
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